Meissner response of anisotropic superconductors 
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The response field of a half-space anisotropic superconductor is evaluated for an arbitrary weak 
external field source. Example sources of a point magnetic moment and a circular current are 
considered in detail. For the penetration depth A<I with L being any other relevant distance (the 
source size, or the distance between the source and the superconductor), the major contribution to 
the response is the A independent field of the source image. It is shown that the absolute value of A 
cannot be extracted from the response field with a better accuracy than that for the source position. 
Similar problems are considered for thin films. 
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PACS numbers: 74.20.-z.74.78.-w, 74. 78.Fk, 74.81.- 



I. INTRODUCTION 

A new experimental technique, the scanning SQUID 
microscopy (SSM), has recently been developed for mea- 
suring magnetic field distributions due to vortices exit- 
ing superconducting samples. 1 Knowing the distributions 
one can, in principle, extract the London penetration 
depth A (either isotropic or anisotropic) and its tem- 
perature dependence. 2 In other implementation of this 
method, one can measure the Meissner response of a su- 
perconductor to a weak external source of static magnetic 
field. Again, the response may provide information about 
the A temperature dependence and its anisotropy. Also, 
one can use tunable external field sources to study ele- 
mentary forces acting upon vortices exiting the surface. 3 
Similar problems are encountered in the magnetic force 
microscopy applied to surfaces of anisotropic supercon- 
ductors. 

There are a few publications addressing these problems 
for isotropic superconductors. 4-6 For anisotropic materi- 
als, however, the response field is asymmetric even when 
the source has certain symmetries, and one cannot use 
methods developed for isotropic materials. To deal with 
the problem, one can utilize the two-dimensional (2D) 
Fourier transform with respect to coordinates x, y of the 
interface, provided the equations for the field distribu- 
tions inside and outside the superconducting half-space 
are linear. This is the case if one adopts the London 
description for the field inside the material. Then, one 
solves the remaining system of ordinary differential equa- 
tions in the variable z, normal to the interface. This ap- 
proach has been developed in Ref. 8 for the problem of 
vortices crossing the superconductor surface. 

Let us consider a source with known field distribution 
h s in the absence of superconductor. In the presence of 
the superconductor occupying the half-space z < 0, the 
total field in vacuum z > can be written as 



h = h s + h r 



(1) 



where h r is the response field which satisfies divh r — 
curl h r = in vacuum. One can look for this field as 
Vtp r with the potential ip r obeying the Laplace equation 
and the zero boundary condition far from the surface. 



The general form of such a potential is 

d 2 k 



tp r (r,z) 



tp r (k)e 



ih-r — kz 



(2) 



Here, r — (x, y) and z is directed normal to the super- 
conducting flat surface at z = 0; ip r (k)e~ kz is the 2D 
Fourier transform with respect to variables x, y at any 
fixed z > 0. The potential (2) is defined only in the up- 
per half-space; hence, there is no problem of uniqueness 
which is in general associated with the description of the 
magnetic field by a potential. 

The field inside the superconductor satisfies London 
equations which read in general anisotropic case as 



hi — A mik eut e-kni h t . n s — , 



x,y, 



(3) 



Here eiki is the unit antisymmetric tensor, my is the mass 
tensor, h t ^ n s abbreviates d 2 h t /dx n dx s . The average pen- 
etration depth A = (AaAftAc) 1 / 3 is related to the actual 
penetration depth for the currents, e.g., along the crystal 
direction a: X a — \^/m a . The masses are normalized so 
that m TO(,m c = 1. 

The problem, therefore, is to match the solutions for 
the field inside and outside the superconductor with 
boundary conditions of the field continuity at the inter- 
face. 



II. RESPONSE FIELD 

Usually in situations of interest, the sample surface is 
normal to one of the principal crystal directions. We call 
this direction c and choose the frame x, y, z as coinciding 
with a, b, c. In this situation, the mass tensor is diagonal 
(m xx = m a , m y y = m b , m zz = m c ), and Eqs. (3) 
reduce to 



hy -\- \ c (H X y hy X ^ X A a (/l 2 ,y 



0. 



),z=0, (4) 



+ ^a{^y,z ~ hz,y),y + {h>x,z ~ hz,x),x — . 

One can replace any of one of these equations by divft, = 
h x , x + h VtV + h ZjZ = 0. It is convenient to replace the 
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third one, and exclude h. z from the first two: 



Solving the system (8) we obtain: 



h X ~ X\\l X ^ XX — ^hx.yy — X\\l X . ZZ + (A^ \\)hy tX y — 0, 

^a,hy,yy ~ X a h yzz + (A c — X a )h XtXy = 0, 

+ hy,y) . (5) 



2„2 



hy X c fl yxx 



Note that the first two equations are decoupled from the 
third. 

We now apply the x,y Fourier transform to Eqs. (5): 

(1 + \\kl + X 2 c k 2 y )h X - + (Ag - \l)k X kyhy = , 

(\ 2 a - Xl)k X kyh X + (1 + \\k\ + \lk 2 y )hy - X^hy = , 

i{k x h x + kyhy) + h' z = . (6) 

Here, hi are functions of k Xl k y and z, and the prime de- 
notes the derivative with respect to z. Hence, we are left 
with the linear system of ordinary 2nd order differential 
equations with respect to the variable z for hi(k, z). 

The solutions are linear combinations of simple expo- 
nentials: 



h l {k,z) = Y J H\ n \k)e q " 



(7) 



The parameters q n and their number are to be deter- 
mined. Substituting each term of Eq. (7) in the system 

(6) we obtain a linear homogeneous system for H^: 



H x (1 + X 2 k 2 x + X\k 2 - X 2 q 



H y i^l 



\ 2 c )k X ky = 0, 



H x (X 2 a - \ 2 c )k x k y + H y (l + \\kl + \ 2 a k 2 y - X 2 a q 2 ) = 0, 
i (k x H x + k y H y ) + qH z = . (8) 

for each n; the superscript n is omitted for brevity. As 
has been mentioned, the first two equations here are de- 
coupled from the third; they have a nonzero solution pro- 
vided their determinant is zero. This gives a quadratic 
equation for q 2 which is readily solved: 



<?1,2 



P 

Q 



P±VQ 
2XlX 2 b 



(9) 



X 2 a + X 2 + X 2 (\lkl + X 2 k 2 ) + X 2 a X 2 k 2 , 
P 2 -4A^(1 + A^ + A^)(1 + A^ 2 ). 



Note that Q < P 2 and both q\ and q 2 are positive; there- 
fore, there are only two positive q's (i.e., n = 1,2) which 
satisfy the requirement of vanishing fields at z = — oo. 

The quantities q determine how the field attenuates in 
the superconductor. We now have to find the "ampli- 
tudes" Hi from Eqs. (8) for each q. It is worth noting 
that solving the homogeneous system of linear equations 
implies, in fact, expressing some unknowns in terms of 
others. To this end, one has to determine the rank of the 
matrix of coefficients for the system (8), choose a proper 
subsystem to solve, etc. The actual procedure might dif- 
fer depending on the situation in question. 



d-X 2 a q 



d = 



k x q(X 2 a -X 2 ) 
1 + Xlk 2 + X 2 k 2 



H, , H, 



d-X 2 q 2 



kyl (XI 



H z 



x\) 

(10) 



for each q of Eq. (9). 

Let us turn now to the field h s of the source. As a 
consequence of div h s = curl h s = out of the source, the 
2D Fourier components of h s are not independent. As 
with the response field, we can look for this field in the 
from h s = V ip s such that V 2 ip s — 0. In our situation, 
the source is situated in the upper half-space z > 0, and 
we are interested in the field h s "under" the source. The 
general solution of the Laplace equation which vanishes 



as z 



-oo is: 



ifi s {r,z) 



d 2 k 



ip s {k)e 



\k-r+kz 



(ii) 



The 2D Fourier components of the source field at the 
interface then read: 

h s a =ik a y s {k) (a = x,y), h s z = k<p s (k). (12) 

Hence, the boundary conditions take the form: 

ik a { V s + ^) = hW+hW (a = x,y), 

k{^ s - V r ) - HP+HP. (13) 

Since the components H a are expressed in terms of H z 's, 
we can solve the system (13) to obtain: 



fc(gi + ff2 - k) - gig 2 (l - l/d) 
k( qi +q 2 + k) + qi q 2 (l~l/d) 

(1) _ k[ip s (q 2 -k)-ip r (q 2 + k)} 

92 - qi 

(2) k[ip s {k-q 1 )+ip r {q l + k)} 

Z 

q 2 - qi 



(14) 



Thus, the response outside and inside the superconductor 
is expressed in terms of the source field (p s . It is worth 
noting that since <fi s (k) can be replaced with h s z (k)/k 7 
the response field can be expressed in terms of the z 
component of the source field at the interface. 

One can see, in particular, that the total flux in z di- 
rection "reflected" by the superconductor is equal and 
opposite in sign to the incident flux of the source cross- 
ing the interface. Indeed, this flux is 

K\ k =o = -kLp r \ k =o = -kip s \k=o = -h s z \k=o ■ (15) 

It is seen from the general formulas (10) that the case 
A a = A;, is singular. The formal reason for this is that 
the first two equations of the system (8) (which we have 
used to express H x and H y in terms of H z ) are no longer 
independent. This situation should be treated separately. 
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A. Small A 

In many situations, the penetration depths are small 
relative to other relevant lengths in the problem such 
as the distance Zq between the source and the interface. 
The characteristic fc's then satisfy k\ <C 1. The relations 
between the response and the source fields then simplify. 
In this approximation, Eqs. (9) give qi = 1/A and q 2 = 
l/X b for A Q < X b . Then, Eqs. (10) yield d = 1 and 

HP=0, H^=-ik v \ a H^\ (16) 
H^=0, H^ = -ik x X b H^. (17) 

Finally, the boundary conditions give: 



y 



1 



X b k x + X a k 2 



fH^ = 2ik yl p° (l - 



), (18) 
X b kl + A a fc 2 - 



k 



*).(19) 



Note that in this approximation A c does not enter the 
result; in other words, the currents along z are small and 
can be disregarded. The results (16)-(19) can also be 
obtained directly starting with the London Eqs. (4) and 
taking advantage of d/dz ^> d/dx a . 

It is also worth observing that in zero order in kX <C 1 , 
cp r = <p s . This means that the x and y components of 
the response field are the same as those for the source, 
whereas h z = —h s z , see Eqs. (2) and (11). In other 
words, in this approximation the response field is the 
mirror image of the source field. 



B. Isotropic materials 

It is readily seen that in this case 

qi=q 2 = VX- 2 + k 2 . (20) 

The first two of Eqs. (8) turn identities, whereas the 
third gives one of H^s in terms of two others, e.g., H z = 
(k x H x + k y Hy)/iq. The boundary conditions then give: 

C. A a = Ab <C A c 



obtain two identical results H x ^k y — Hy^k x — which 
together with the third equation yield: 

H li ) = iQ^ H li )j H W = Wh H W. (23) 



Doing the same for q 2 , we obtain from the first two equa- 

(2) (2) 

tions H x ' k x + H y ' k y — which is compatible with the 
third only if Hi 2) = 0. Hence, we have: 

(24) 



Jj(2) =0 H (2) = _^ H (2) 

rvy 



The boundary conditions (13) yield: 



qi+k 



qi+k 



Note that for this case, q 2 along with A c drop off the 
result for any source. In other words, the response is as 
if the superconductor were isotropic with the penetration 
depth X ab . 



D. A c = A a < Aj, 

This is the case of the screening by the "side surface" of 
a uniaxial crystal like BSCCO. Our notation is, however, 
differs from the commonly used (for standard notation, 
we should have replaced in our formulas A c and A a with 
X ab and Afc — ► A c ). We then obtain using Eq. (9): 



qi = VXa 2 + k 2 , <?2 = Y 1 + X\k 2 + X 2 k 2 J X b . (26) 
With q = qi, Eq. (10) gives 

HW=i±H?), HW=i±+£&HW t (27) 
qi y K k y<H 

and for q = q 2 



k 2 ' 

r« »„a rr(2) 



(28) 



Finally, Eqs. (14) give tp r , Hi , and Hi ' in terms of the 
source field ip s . 



This is the case, e.g., of BSCCO crystal with the ab 
plane being the surface. Equations (9) yield: 

qi = yjK? + v, 92 = v/A- 6 2 + 7 2 fc 2 , (22) 

where 7 = X c /X ab is the anisotropy parameter. Substi- 
tuting qi in the first two equations of the system (8), we 



III. EXAMPLES OF SOURCES 

To apply formulas derived above for the response fields 
outside and inside superconductor, one needs the Fourier 
transform (p s (k) of the source field at the interface. Below 
we provide examples for which ip s (k) can be calculated 
analytically. 
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A. Point magnetic moment 

Consider a magnetic moment fi situated at the hight zo 
above the superconductor. The corresponding potential 
(in the absence of the superconductor) at z = is 



where fi± is the in-plane part of the magnetic moment, 
and (3 is the angle between and x. For A a = A^, the 
position of (i± in the xy plane is arbitrary, still there is 
a torque which tends to rotate fi out of the z direction 
and to place it in the xy plane. 



H R [i z z - fi ■ r 



R 3 



(r 2 + z 2 ) 3 / 2 



(29) 



B. The source as a current loop 



Here, R = (x,y,z — z ) is the radius-vector originating 
at the source; the first minus sign is due to the definition 
h s = \7tp s . The 2D Fourier transform is: 9 



(30) 



The response field can now be calculated with the help 
of Eqs. (14). In general, this can be done numerically; 
for A <C zq, the analytic evaluation is possible. 

As an example take the moment directed along z 
above the flat isotropic superconducting surface; for the 
isotropic thin film this problem has been considered in 
Ref. 5. According to Eq. (18) for the isotropic case 
tp r = ip s (l — 2kX). Transforming back to real space we 
obtain: 



(p r {r,z) = 



2A- 



2Z 2 



R 



\f~r 



Z 2 , 



R 5 
Z 



Zq . 



(31) 



Here, the first term is the field of a moment — /x at z = 
—zq, i.e., of the image source. The second term is the 
field of the proportional to A magnetic quadrupole at the 
same point. 

For the magnetic force microscopy, the quantity of in- 
terest is the interaction energy which is given by 



£ = - (J,- h r (0,0, z ) 
d 2 k 



(32) 



(27T) 



Substituting here Eqs. (18) and (30) and integrating we 
obtain: 



£ = 



+ 



Az 3 



32z\ 



1 + 



3(A a + A b ) 



4z 



+ 



t4 + 14 

8z$ 



[nl(3X b + X a ) + fi 2 y (3X a + X b )]. (33) 



In addition to a repulsive force —d£/dzo, the magnetic 
moment fi experiences a torque, because the energy de- 
pends on the moment orientation. It is readily seen that 
if the position and the value of the magnetic moment are 
fixed, the minimum of £ corresponds to situated in the 
plane xy and parallel to the direction of largest A. The 
z component of the torque is easily evaluated: 



3Mi 
16z 4 



(A a - X b ) sin 2/3, 



(34) 



Let the source be a circular current of a radius a situ- 
ated in the plane z — z - If a and zq are of the same order 
of magnitude (practically, they are both of a few-micron 
size), modelling of the source by a point-size magnetic 
moment does not suffice. The scalar potential of the field 
created by a loop in the plane z = reads: 7 



f' '(r) 



dS-R 

R? 



(35) 



The source current I flows in the counterclockwise di- 
rection relative to the z axis so that an area element 
dS = dS z. R is the radius- vector from this element 
to the interface point (f, 0), and the integral is over the 
area of the current contour. The position of the ele- 
ment dS = d 2 r'z in our situation is (r',z ) so that 
R = r — r' — z z, and 



^ s (r,0) = 



Iz 



d 2 r' 



[(r - r') 2 + Z 2]3/2 ' 



(36) 



For the circular loop of the radius a, the integral here is 
over the circle area. Comparing this with Eq. (29) we 
see that the source field can be considered as created by 
magnetic moments distributed uniformly over the loop 
area with the density Iz/c so that the total moment of 
the loop is wa 2 Iz/c. 

The 2D Fourier transform of z [(r - r') 2 + z 2 ]~ 3 / 2 with 
respect to r is 2ire~ kz ° e ~ lk r , see Eq. (30); therefore, 



V s (fe) 



- 2i:e 
c 

4w 2 Ia 

ck 



■kz 



d 2 r'e 



I 

fcZ0 Ji(fca 



-ik-r' 



(37) 



According to (25) the 2D Fourier transform of the re- 
sponse potential for isotropic superconductor is: 



<p r (k) = 



4ir 2 Ia q - k 



-kz 



Ji(ka) . 



ck q + k 

The z-component of the response field follows: 



(38) 



h r z (r,z) = 
2ixla 



'i!' 



v - 1 



f 

Jo 



q + k 



e -H*+*o) j 1 (ka)e ih - r 



dk k 



q — k 
q + k' 



-k(z+z ) 



Ji{ka)J {kr) . (39) 



Further, one can evaluate the response flux through a 
flat probe placed above the superconductor. If the probe 



■5 



is a circular loop of a radius a p with the center at r = 
at the height z p < z 0l the flux is given by 



4n 2 Iaa 



~ [°° dkS —T e ~ Kz ° +Zp)j i(ka)Ji(ka p ). 
Jo q + k 

(40) 

All formulas for the isotropic case have been worked out 
earlier by Clem and Coffey making use of the cylindrical 
symmetry of the problem. 4 



C. Interaction with vortices 

It is of interest to evaluate the force acting on the vor- 
tex tip by the screening currents in the sample created 
by a circular current as the field source; experiments for 
which this is relevant are described in Ref. 3. A similar 
problem for isotropic films and the magnetic moment as 
a source has been considered in Ref. 5. We do this for 
materials isotropic in the xy plane, for which the force 
depends only on the distance r from the loop center along 
with the loop hight z . One can calculate F x (x, 0; z ) and 
in the result replace x with r: 

4>0 f° 4>o f° 

F x (x,0) = — dzj y = — I dz(h x . z - h z , x ) 

C J — CO J CO 



With the help of Eqs. (22), (23), and (37) we obtain 
<j) Ia f°° dk k e- kza 



F r (r) 



cKb Jo 9i (9i +k) 



Ji(fca)Ji(fcr) . (42) 



Using 1/A 2 = q\ — k 2 and the substitution t = ka, we 
write: 



F r . = - 



2ca 



f°° /A \ 
J dtG{-t)te-* t J 1 (t)J 1 (r't), 



A 2 ,\-V2 



G = l-±t(l + ^A 
a V or / 



(43) 



where z' = z /a and r' = r/a. Usually, the parameter 
A/a <C 1. Besides, due to the factor e" z t Ji(t), the re- 
gion contributing to the integral is < t < min(l/z', 1) 
because of the oscillating Ji(t) at large t (unless r' w 1). 
Then, one can expand G in powers of At/a and keep only 
the linear term: 



f- CO 

fo = / dtte-*' t J 1 {t)J 1 {r>t), 
Jo 

/•OO 

h = / dtt 2 e-*' t J 1 (t)J 1 (r , t) ; 
Jo 



(44) 



The integrals here can be expressed in terms of the hy- 
pergeometric functions convenient for the numerical eval- 
uation, see Appendix A. 



One can also define a potential energy U(r) such that 
F r = -dU/dr: 



TT ( X \ 



u 



dte- z,t Ji{t)J {r't) , 



(45) 



poo 

= / dtte- z,t Ji(t)Jo(r't), 
Jo 



This energy is defined so that U(oo) = 0; its value at 
the origin is 



U(0) = 



2c 



1 - 



Xa 2 



{z 2 + a 2 y/ 2 (z 2 + a 2 )3/2 



(46) 



Thus, the source loop creates a potential well of the depth 
|[/(0)| or a barrier of the hight |£7(0)| for vortices under- 
neath depending on the current and vortex directions. 
This opens an interesting possibility for studying the be- 
havior of vortices (or antivortices) in tunable potentials. 3 
In a similar manner, one can evaluate interactions of vor- 
tices with other types of sources. 



D. Accuracy of the SSM determination of A 

Magnetic fluxes of the response field, in particular, 
$ 2 = J d 2 r h r z (r, z p ) (the integral is over the area of a 
pick-up coil placed at the hight z p above the supercon- 
ducting surface) can be measured with high accuracy for 
a given source, given geometry of the coil, and a known 
height z p . In principle, this leads to a possibility to mea- 
sure the penetration depth. However, the accuracy of 
this determination in bound by the accuracy with which 
the heights zq and z p are known. To demonstrate this 
consider the response field of an isotropic material for 
small A's: 



h r z (r,z p ) 



d 2 k 
W) 2 



r ( i,\Ak-r—kz v 



k(p r (k)e 



= - I kp s (fc) (1 + 2\ky k - r - kz * . (47) 

J (2tt) z 

If z p varies by Sz p , the response field variation is 

sk = sz p J -0^ fc V(fc)(i + 2\ky k r - . (48) 

If only A varies, we have: 

8h r z = -2SX J -^^k 2 if s (k)(l + 2Xk)e lk - r - kz ^ (49) 
in other words, 

(50) 



6K „ 5h r , 
SX 5z p 
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Therefore, the accuracy of extracting A from the data 
on h z cannot be much better than the knowledge of z p 
(the same is true about the source position z n ). The 
latter is usually known within a fraction of a micron. This 
is a severe restriction upon the accuracy of the absolute 
determination of A. Still, in principle, SSM allows one 
to determine accurately the temperature dependence of 
A (for fixed z p and z n ). 



IV. THIN FILMS 

The Meissner response of superconducting thin films 
can be probed by the SSM method in yet greater detail 
than that of bulk samples. The films often have a large 
Pearl length A = 2\ 2 /d (d is the film thickness) exceeding 
substantially the size of the SSM sensing loop, making 
the SSM measurement to a local probe. Formally, the 
problem of a thin film in a field of an external source 
is simpler than that of the superconducting half-space, 
because in the film case there is no "internal problem" 
to solve; instead, the film provides a boundary condition 
for the outside field distribution. 

Consider a film in the x,y plane made of a uniaxial 
material with the c axis at an angle 9 to the film normal 
z. We write the London equation (3) for i = z, h z — 
4ir\ 2 (m x iji tV — m y iji tX )/c = 0, and integrate it over the 
film thickness: 



2tt 

h z - —A(m xx g Xiy - m a g y ^ x ) = , A 



2A 2 
d 



(51) 



where g is the sheet current. Note that only two com- 
ponents of the mass tensor, m xx — m a cos 2 9 + m c sin 2 9 
and rriyy — m a , determine the film anisotropy. 

Using the relation of sheet currents g to the tangential 
components of the response field, 



2tt 



-^(+0), ^g y = K(+0), (52) 



(+0 denotes the upper film face; the tangential compo- 
nents satisfy /i£(+0) = — 0)), we obtain for z = +0: 

h s z + h r x +A(m xx h r yty + m a h r x>x )=0, (53) 

for more detail, see Ref.8. Further, since h s z +h r z = k((p s — 
ip r ) we obtain: 



<P' = 



k<p s 



k + A(m xx k 2 + m a kl) 
In particular, for 9 = we have 



<P' = 



z = +0. (54) 



2A 



1 + A afc < K=m a K=^. (55) 



This result holds also for isotropic materials where 
m a = 1. In this case, one can readily obtain the field 
h z for the circular current (37) at the height zq above 



the film; its 2D Fourier transform for < z < z is given 

by 



h r z (k,z) = -k(p r e' kz 

At: 2 la J^ka) ^ k{z+Za) 
c 1 + kA 



(56) 



This field can be measured by SSM. 10 

For completeness, we write down also the field on the 
opposite film side, i.e for z < 0, where the response po- 
tential is given by 



(p r (r,z < 0) = 



d 2 k 

(2n) 2 



r (h\ „ik-r+kz 



if r (k)e 



(57) 



The London boundary condition (53) should now be writ- 
ten in terms of the field components at z = — 0: 

h s z + h r z - A(m xx h r y y + m a h r x x ) = , (58) 

which yields the 2D Fourier transform of the response 
potential (54) with the minus sign. Proceeding as above, 
we obtain for the isotropic case: 



h r Mz<0) = _^£m e H*-*o) 



(59) 



The force acting upon a Pearl vortex situated in the 
film at a radial distance r from the current ring center is 
readily evaluated: 

F r {r)=-^— — -Ji(fca)Ji(fcr)dfc. (60) 



c J 1 + kA 
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APPENDIX A 

The integral u of Eq. (45) is known (see Ref.9, 
6.612.3): 

poo 

I(r,z) = / e- zt Ji(t)Ji(rt)dt 
Jo 

= 25/2^3/2 2Fl (l'4 ;2; ^)' (A1) 



7 



where the primes by z and r are omitted, and u = [z + 



l)/2r. Therefore: 



f 

Jo 



te~ zt Ji(t)Ji(rt) dt = - 
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(A2) 



3zr -3/2 



2 7/2 u 5/2 



^(l'4 ;2; ^) + 8^ 2Fl (4'4 ;3; ^) 
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